(13) 



where R = rank(Bo). Finally, a balanced realization of the system under question is obtained 
by pseudo-inverting various submatrix components (block 108 of FIGURE 1) as follows: 



Since R « L, the above model represents a reduced-order realization of the original 
system. Note that the realization is optimal in that it is balanced between inputs and outputs. 
10 However, the total number of samples taken is N ( x (M+ 1), which increases proportional to 
the number of inputs. Also, for an accurate system identification Ho must have sufficient 
columns and rows, i.e., N f x (M- 1), >R and N 0 >R. 

For a very large data set Y" with many time steps and a large number of inputs, the 
Eigensystem Realization Algorithm/Data Correlations (ERA/DC) method may be preferred. 
15 As described, for example, in Juang, J.N., Applied System Identification. Prentice Hall 
Englewood Cliffs, New Jersey, 1994, the ERA/DC method may be used to compress the 
amount of data and reduce the computation time required for the SVD of the Hankel matrix. 

Although desirable results have been achieved using the prior art computational 
methods, there is room for improvement. For example, if the unsteady CFD model is driven 



5 



D = Y° (N 0 x 1) 

C ~ U*£}/ 2 (N 0 xR) 

B ~ the first N t columns of E^ /2 V£ (R x iV<) 

A ~ S^U^VkS^ 2 (RxR) 



(14) 
(15) 
(16) 
(17) 
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With reference to FIGURE 2, the SCI/ERA method 200 proceeds as follows. First, at 
a block 202, individual pulse responses are sampled for the first two time steps: 

Y° = far? y 2 ° • ■ • yfc] (18) 
Y 1 = [y! y\... y y (19) 

Next, at a block 204, an SCI is constructed as follows: 



Ni 

h sci = E b i r " (for states) (20) 
t=i 

Ni 

d n SCI = Edir? (for outputs) (21) 

where 

r t n = a sequence of arbitrary numbers (22) 



10 To ensure independency of the inputs, it is desirable that the signals be as 

uncorrelated as possible. In an ideal case they would be statistically uncorrelated random 
signals, i.e., Q/m) =£/r " r\ m J = 0 for i but they are hard to construct for numerical 
analysis. 

Subject to the SCI, at a block 206, the SCI/ERA method 200 samples the system 
15 response y" for n = 0, 1, 2, . . . , M t to obtain: 

y c n o = cx" 

= y n -Ey?^ n ( 23 ) 
= y n+1 - £ yf r? +1 - £ y] r? (24) 

t=l i=l 



i 
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